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1.   Introduction 

The  deforraatlons  of  an  elastic  cylinder  are  deterrained 
if  a  solution  of  the  equations  of  the  exact  tliree  dimensional 
linear  theory  of  elasticity,  hereaftei'  referred  to  as  the 
exact  theory,  can  be  obtained  vjhich  satisfies  the  prescribed 
force  and  displaceraent  co.iditions  on  the  surfaces  and  edges 
of  the  cylinder.   Since  it  is  di.i' f lo.lt,  in  general,    to  con- 
struct this  solution  aporoxiaiato  theories  have  been  proposed 
(for  cxaraole,  see  ref.  1-5).   These  theories,  called  thin 
shell  theories,  are  derived  frova  setc  of  a  priori  assuraptions 
concerning  the  dir^tribution  of  stresses  in  the  shell  and  con- 
sist of  a  determined  system  of  differential  equctions,  boundary 
conditions  and  stress  dis-^lacement  relations.   These  asr;uiaptions, 
x-'hich  are  often  mutually  contradictorjr  and/or  inconsistent  viith 
the  exact  theory,  can  be  accurate  only  if  the  shell  is  suf- 
ficiently thin.   In  spite  of  these  inadequacies  reasonable 
aporoxlm.ations  are  obtained  in  la&nj   crses.   liowever,  it  is 
difficult  to  ascertain  the  relative  accuracies  of  these  thin 
shell  theories,  x-ihich  vary  in  complexity,  and  to  establish 
their  relation  to  the  exact  theory.   It  is  not  our  Intention 
to  oresent  a  complete  or  critical  rev^ev  of  the  e::lsting 
literature.   For  discus^;,  ions  of  tliis  nature  the  reader  is 
directed  to  references  2,  5,  6,  7,  8  and  references  contained 
therein. 


It  is  reasonably  conjectured  that  the  thin  shell  theories 
rcay  -"-ield  Inaccurabe  results  in  the  neighborhood  of  boundaries 
(for  example,  nee.r   the  edges  oj.'  the  cylinder  or  in  a  region 
adjacent  to  a  hole  in  the  cylinder  -all)  and  for  "thicker'' 
shells.   Attempts  have  been  made,  by  employins  "less  restrictive" 
a  priori  assujnptions ,  to  obtain  appropriate  thick  shell  theories. 
Some  of  these  theories'"'  '        are,  in  a  sense,  generalizations 
of  Reissner's  theory  5or  thick  plates.     See  ref.  8  lor  refer- 
ences and  a  discussion  of  so^i'e  recent  Kork.   However,  the 
nature  of  the  assu^mtions  a';^ain  nakes  it  difficult  to  estimate 
the  accuracy  of  such  theories  and  their  relation  to  the  e:cact 
theory.  •  ' 

It  is  the  pur-pose  of  this  orper  to  rectify,  to  sotie  e:.tent, 
the  inadequacies  in  these  approachec-to  shell  theory  by  system- 
atically deriving  fror.i  the  exact  theory  approximate  thin  and. 
thick  shell  theories.   This  includes  boundary  conditions  as 
well  as  differential  equacioriS  and  stress  displacement  relations. 
VJe  only  consider  the  rotationally  symmetric  deforrjations  of  a 
cylindrical  shell  of  radius  R  and  uniform  T'all  thickness  2h, 
However  the  procedures  that  are  eirroloyed  may  be  suitably  ex- 
tended to  unsymmetric  defornations  and  to  shells  of  otl:cr  shapes. 

The  method  of  derivation  begins  by  assuming,  after  aporo- 
nr lately  scaling  the  independent  and  dependent  variables,  '  that 
each  stress  and  disolacement  cOi,iponent  can  be  represented  by 
a  nower  series  expansion  in  s,  a  small  parameter  which  is  a 


'   The  proper  scaling  is  a  crucial  point  in  the  method.   It 
deter-iines  the  type  and  "width''  of  boundary  layers  under 
cons  ider at  ion . 


positive  povjer  of  h/Pi  •   These  series  are  then  substituted  into 
the  exact  equations  of  elasticity  and  the  exact  boundary  con- 
ditions on  the  inner  and  outer  sur'faces  of  the  shell.   By 
equatinfc  coefficients  of  similar  poi/ers  of  e  ir  these  equations 
a  sequence  of  algebrsic  and  differential  equations  is  obtained 
to  deter-Jiine  the  coefficients  in  the  expansions  of  the  ctresses 
and  displacetiients .   The  first  of  this  set  of  differential  equa- 
tions and  stress  displacement  relations  are  identical  with 
those  of  the  classical  thin  shell  theory  of  Love,"^   .Higher 
order  terras  yield  'thick  shell'  covrectiors.   I.'owever,  these 
expansions  do  not  satisfy  the  boundaiy  conditions  along  the 
edges  of  the  shell  nor  can  they  represent  the  stresses  and 
displacements  in  a  region,  called  the  boundary  layer,  adjacent 
to  the  edge.   To  deter-jiiine,  fro:n  the  exact  theory,  appropriate 
boundary  conditions  for  the  differentlDl  equations  and  to 
exaraine  the  stresses  in  the  neighborhood  of  the  edge  vje  eraploj?- 

a  boundary  layer  expansion  technique  related  to  the  one  pro- 

1?  l'^ 

posed  by  Priedrichs  "  and  later  Priedrichs  and  Dressier  "  m 

a  study  of  t^xe  bending  of  nlates.   The  expansion  and  scaling 
procedures  thPt  we  employ  differ  somewhat  fron  those  in  refs. 
12  and  13.   in  so:'ne  sense  oux'  uiethod  is  a  generalization  of 
the  one  given  in  these  references.   The  stresses  in  the  neigh- 
borhood of  the  edge  are  then  systeidatically  approxiniated  by 
eraplovlng  the  solutions  of  a  sequence  of  boundary  value  prob- 
leras  for  the  bihariaonic  equation  in  a  serii-lni  inite  strip. 
VJe  then  show,  in  terms  of  the  definition  given  In  Section  3, 


that  there  is  only  one  scaling  of  the  independent  variables, 
i.e.  one  boundary  layer,  that  leads  to  a  "reasonable  shell 
theory" . 

A  portion  of  our  results  are  obtained  in  a  so^iiewhat  less 
systematic  manner,  in  an  unpublished  paper  by  Johnson  and 
Reissner,  "^  hox^jever  the  boundary  layer  is  not  analyzed  in 
this  reference;   so  that  the  boundary  conditions  for  the 
differential  equations  are  essentially  assuiaed  rather  than 
systematically  derived  from  the  exact  theory,  and  no  procedure 
is  given  to  approximate  the  stresses  near  the  edge.   The 
relation  of  the  results  of  ref .  l5  to  those  of  this  paper 
will  be  made  clearer  in  the  following. 


2.   Formulation 

To  study  the  deformations  of  a  circular  cylindrical  shell 
it  is  convenient  to  introduce  a  cylindrical  coordinate  systeiii 
(r,e,x).'   VJe  consider  a  cylindrical  shell  of  radius  a,  thick- 
ness 2h,  and  length  L  as  a  three  dimensional  body  bounded  by 
the  coaxial  cylindrical  surfaces  r  =  R  +  h  (called  the  outer 
surface)  and  r  =  R  -  h  (called  the  inner  surface),  and  the 
parallel  planes  x  =  0  and  x  =  L.   The  annular  regions,  x  =  0,L, 
R-h<r<R+h,  are  called  the  edges  of  the  shell*   The 


It  is  customary  in  shell  theorj'-  to  introduce  the  symbol  x 
instead  of  z  to  denote  the  axial  coordinate.   The  sjnnbol  z  is 
usually  reserved  for  representing  distances  along  the  normal 
to  the  middle  surface  of  the  cylinder. 


surface  r  =  R  is  referred  to  as  the  middle  surface.   V'e  should, 
for  consistence^  i/ith  the  scaling  orocedure  introduced  in  eqs. 
(1),  place  the  origin  of  coordinates  at  the  center  of  the 
c-'linder  vdth  tlie  ed^es  then  given  by,  x  =  +  L/2,   Iio^^Jever, 
it  is  convenient  for  our  purposes,  and  there  is  no  loss  of 
generality,  to  hrve  the  orir/in  coincide  i-.'ith  the  end  of  the 
shell. 

It  is  aseunied  that  the  cylinder  is  constructed  from  a 
homogeneous  and  isotropic  elastic  :aaterial  and  that  the  deform- 
ations are  sufficiently  snail  so  that  the  linear  theory  of 
elasticity  is  valid.   The  external  x'orces  acting  on  the  cylinder 
are  nor  'al  forces  on  the  inner  and  outer  surfaces  and  normal 
and  tangential  forces  applied  to  the  edges.   All  forces  are 
distributed  in  a  rotationally  symmetric  manner.   If  it  is 
assumed  that  the  resulting  deformations  of  the  cylinder  are 
also  rotationally  syi.rmetric,  then  the  only  non-vanishing  dis- 
placements, [X   and  O  in  the  x  and  r  directions  respectively, 

and  the  onl^^  non-vanishing  stresses,  -t      ,  ^r-r^,  '^"      ,    and  'c      , 

o         *   rr'   r@'   XX'      rx* 

are  functions  of  x  and  r.   I'e  define  the  parameter,  e,  as. 


1 


=  (&)^  ,  (1) 


where  b  ie  a  positive  integer  and  introduce  the  dimensionless 
variables ; 


r  -  ^^  -   X  „  _  r-R  _  r-R  ,  ^^  x 

^-j— S7T5    '  ^-~b---E-    '  (2a) 

R(|) 


cy{Z,K}t)    =  -^ ,        -7^(5, ^,'S)  =  --^1. ,  I   (2b) 


-  (r.,^;e)  =  ^ ,      '^-,^(?.^;e)  =  J^ / 


Here  E  Is  Young's  modulus  end  a  is  a  non-necative  Integer 
such  that, ' 


X  ^  b  -  a  >  0  .  (3) 

VJith  eq,  (3)  the  possible  (a,b)  coribinations  are  li:iited  to 
the  sector,  tt/Ij.  <  a  <  n/Z,    in  the  a,b  plane  where  a  is  the 
polar  angle. 

The  independent  variables,  <4  and  ^,  are  respectively  the 
scaled  axial  coordinate  and  the  rcaled  "thickness"  coordinate. 
Prora  (2a)  the  equations  for  the  outer  and  inner  surfaces  are 
given  respectively  by,''  ^  -■   +1,  ^  =  -1.   An  interpretation  of 
the  scaling  procedure  is  given  in  Section  Ij . 

In  terms  of  the  diinensionlsss  variables  (2)  the  equations 
of  the  exact  theory  with  rotational  symraetry  are  written  as  ; 


'   In  the  follovjing  section  it  is  shown  that  values  of  a  and  b 
are  deternined  by  the  order  of  magnitude  of  the  external  forces. 
The  inequality  in  (3)  is  a  consequence  of  the  scaling  (2a)  and 
the  "stretching"  procedures  introduced  in  Section  I4.. 

-h 

It  is  possible  to  obtain  a  riore  general  scaling  procedure  than 
(1)  and  (2a)  by  replacing  the  exoonent,  l/b,  in  (i)  by  an  arbi- 
trary non-negative  integer.   The  coiralexity  of  the  calculations 
are  then  considerably  increased  and  no  ne\<!   information  is  obtained, 


Equilibrium  Equations, 


Sti-ess-Displaceraeiit   Relatious    (Hooke's   Lavj), 


If  =  e^[o^-  v(T^+,r  )]      ,  w  --=   (1+  e'^?)(<.g-v(<T^+r  )]    , 


V 


(5: 


Compatibility  Equations, 


?\   -2, 


Ao-  +  li;  ^"^'- 


X  '  1+v      ,;2      "    '  \ 

1 

?■ 
/ 
I 


(6 


where    v   is   Polsson's   ratio   and, 

b 
1+e^  z       y       X 


To  coiirolete  the  formulation  of  the  exact  theory  we  obtain, 
from  the  specification  of  the  applied  external  forces  the  fol- 
lowing boundary  conditions: 


0-^(5, -l;e)  =  P3_(^;e)  , 


(7a) 


V(7b) 


R(|)  ri(|) 

The  functions  p^C^je  )  ,p-^  (£>;e)  ,o(^;e  )  ,~^(^;e)  and~2(^j^^  ^^® 
arbitrary  to  the  extent  that  they  form  a  system  of  forces  in 
equilibrium  and, 

Eq.  (8)  im.plies  that  surface  forces  which  only  stretch  or  only 
compress  the  cylinder  v;all  are  not  permitted.   This  restriction 
is  not  essential,  it  merely  sir.iplifles  the  ensuing  calculations. 

A  complete  mathematical  formulation  of  the  elastic  prob- 
lem for  the  rotationally  symmetric  deformations  of  a  circular 


v-j_(+l)  =  -.'2(11)  =  0  for  continuity. 


cylinder  (ior-iiulftion  A)  con-iat^c  of   the  equilibrium  equations 
V-i-) ,    the  streas-dlsplacerient  relations  (5)  and  the  boundary 
conditions  (7).   A  second  and  equivalent  formulation  (form- 
ulation 3)  consists  of  the  equilibrium  equations  (.'4.),  the 
compatibility  equations  (6)  and  the  boundary  conditions  (7). 
In  the  follovilnr^^:  section  x'e  employ  formulation  A;  i/hlle  in 
the  boundary  layer  invest i.^ation  In  Sections  Ij.  and  5  It  is 
convenle:it  to  use  for^iulatlon  3. 

Section  6  contains  a  summary  of  the  shell  tlieory. 


3.   The  interior  Problems 

To  obtain  in  a  systematic  rp.anner  from  the  exact  theory, 
(ij.)  -  (7),  approximate  '"thin"  and  ''thick"  shell  theories  we 
adopt  form.ulation  A  of  the  previous  section  and  ascurae  that 
each  stress  and  dlsplace'sient  component.  Indicated  by  the 
generic  sym.bol  <r(4,c^;  e ) ,  can  be  represented  by  a  power  series 
in  e : 

r{^>K;B)^--  ZI^o   -  (^»^)e^  >  (9) 

where  rr     =  0  if  1  <  0.   The  coefficients,  ■t''"(^5  ?),  are  referred 

to  as  the  interior  stress  coefficients  or  the  interior  dis- 
placement coefficients,  whichever  the  case  may  be.  Ve    assume 
that  the  prescribed  external  forces  in  (7)  can  be  written  as: 


10 

Expansions  of  the  form  (9)  for  the  stresses  and  di;jy)lace- 
ments  are  substituted  with  (10)  into  (/j.),  (5)  and  (7a). 
Equating  of  coefficients  of  the  same  power  of  e  in  each  of 
these  equations  results  in  systerris  of  algebraic  and  differential 
equations  and  surface  boundar:,-  conditions  for  the  interior 
stress  and  displacement  coefficients.   Froi'i  the  coefficients 
of  e^  we  obtain, 

n    ., -.n-A 


^^5  +  _4-^  +  7"'  ■.-^,.  =  0  ,  (lib) 


=  a-.   -  vGt—  +<--)  , 


."  =  0^  -  v(a;;  +o-^)  -H  [.-^-^  -  v(T^-^  +.7-^-^]^   j 


(12a,b) 


__  =  ,r^      -vCt^_  +0-^   )  ,   ■:t^  +  -~-^--  -  2(l+v)o  ^^  ,   (12c, d) 

^-^  -^    {"  ^. 

.4(?,11)  =  0  ,   Scr^(^,l),a^(?,-l)J  =  ^p^(?),p5(£,)>  .   (13a, b) 


Here  the  binonisl  expansion  of  (1+e''^)"'  in  g(^)  has  been  used 


and, 

i+Tj+TTb^n 


''-  'a^  ™  ^  2l   (-1)^*  ^^'a^  ,        i,j  >  0  . 


11 


We  shall  now  obtain  from  (11-13)  appropriate  cliff erei-itial 
equations  and  etresc  displacerasnt  relations  for  the  Interior 
stress  and  displacement  coefficients.   To  begin,  eqs.  (12a, b) 
are  solved  for  J      and  CTq, 

^"^'^    n,.   V   n    V      r^.n-b  ./.Ji-b.  n-b  n  i  ^   i 


n    1  /3u    ,   n\,   V  n    v  r-rn-b  ,, /^Ji-b  ,  n-b -, -i  ^ 
1-v   ^  1-v 


X 


V  (li|) 


1-v     -^  1~V  ,         y 


Then  frora  (3)  and  (12d)  vlth  n  =  a  \.e    see  that  either. 


9?      ^^? 


f  2a  -  b  >  0  ,  (Case  I) 


m5) 


u^  =  U^(?)  if  2a  -  b  <  0  ,  (Case  II)  ,j 


where  we  have  use6,o''^   =  0  if  1  <  0.   Case  II  leads  to  "membrane 

theories''  and  will  not  be  discussed  here,   VJe  will  therefore 

restrict  our  attention  to  Case  I,   The  special  case,  2a-b  =  0 

(Case  I-A),  is  considered  first. 

J. 
Case  I-A;   In  this  ccse  the  coordinate  scallncs,  (2a),  become. 


Rh 


'   Johnson  and  Reissner  ^tnplo.-ed  a  scaling  equivalent  to  (16) 
but  arrived  at  it  throu'^h  a  prio.ri  ■^^byslcal  reasoning.   u'e  I'-iill 
show  later  that  this  is  only  scaling  that  gives  a  reasonable 
shell  theory. 


12 


and,  A  =  b  -  a  =  a.   It  then  follows  from  (lib)  that, 

——   =  0   If  n  <  a  ,        -—■   =  0   if  n  <  Pa  . 


Integrating  these  expressions  and  using  the  surface  boundary 
conditions,  (13),  we  find  that^ 

rr'^     =0   if  n  <  a   ,      q-     =  o^  ~  p""  =0   if  n  <  2a  , 
zz  '       z   '  o    -I  * 

(17a, b) 

vjhere    (6)    has   been  employed.      Ve    conclude   frora    (17b)    and    (ICa) 
that   surface   loads   i.hose   orders    of  maj^^nltude   are   less    than, 
e"^"^,    and   satisfy    (8),  crj.nnot  be    prescribed. 

It   follows    iriiined lately  from   (12c)    that, 

^i'{^,K)   =  ^'^iO  if   n  <   2a      ,  (l8) 

and  from  (12a, c,d)  that,  within  a  rigid  body  motion, 

u^(5,;C)  =  0        if  n  <  a  . 

Integrating  (12d)  and  using  (17^^)  £uid  (18)  there  results, 

u^(i,^)  =  -  ^^^-^  K  +   U^('^)    if  a  <  n  <  3a  .    (19) 

Substituting   this   i-esult,    (17b)    and    (l8)    into   eqs.    (iZ-f.),    the 
first   2a   interior   stress-dlsplaceinent   relations    for   the   coef- 
ficients  of  (T       and  o\   can  be   vjrltten  as, 


13 
a^(5,^)  =  3:^(5  )->-C(^J^  ,  o'ii^^K)   =  S^^(4)+-2(C)^  if  n  <  2a  ,  (20) 

vhere. 


X      i^v^^     '^  "^      1-v      ^ 

^n..^  _    1   dV(C)        en,,.  _  _  _y__  dV^(^) 
^       l-v^   d?"-  ^       1-v^   dg^        y 

a^^   and'T^  are  therefore  knox:n  if  n  <  2a  xjhon  V.'   and  U    are 
deterrf^.lned  for  n  <  2s.   We  novj  obtain  using  (11)  and  (13) 
appropriate  differential  equations  bo  deterinine  these  quantitie: 
as  well  ac  stress  displacement  relations  i'or  O"^^^  and-j'^. 

Substituting  (20)  into  (lib)  with  n  in  the  range, 
a  <  n  <  3a,  integrating  with  respect  to  ^  and  employing  (13a) 
we  obtain  with  the  aid  of  (17a)  s 


'   The  coefficients  defined  in  (21)  arc  related  to  the 
conventional  stress  i^esultants,  II  (x)  and  n^(x),  and  couples, 
M  (x)  and  Ii(_Jx)  of  the  classical  shell  theories.  '     If  each 
of  these  quantities  is  expanded  in  a  power  series  in  e  then  it 
is  easy  to  shov/  that, 

n"(x)  ^      n2(x) 


X 


i^)  =  ~^r   '  s2(U  =-^K-   ' 


3Ii"(x)         _.^     3m"(x) 
^       2Eh  "       2Eh'^ 


Ik 


7. 


If  n  <  2a  , 


(22) 


and, 


'■"x'^-?>  =  I  -c-'^-?"' 


if 


<  n  <  3a  , 


(23, 


Eq.  (22)  yields  one  of  the  differential  equations  to  determine 
ea.ch  of  the  first  2a  interior  stress  coefficients;  uhile  (23) 
gives  the  stress  displacement  relation  for  the  flr^t  2a  non- 
vanishing  r\      .   To  obtain  the  second  differential  equation 

and  stress  displacement  relstions  foi''  o   we  consider  (11a) 

z 

with  n  in  the  range,  2a  <  n  <  [|.a.   Substituting  (20),  (23)  and 
(17b),  integrating  once  viith   respect  to  c,   and  inserting  this 
result  into  (13b)  we  find  that. 


^2^n-2a 

X      oon-2a  „  -, ,  n  ^n^ 

5 3Sg     -  ^(Pj-Pq) 


d£>- 


if  2a  <  n  <  i|a  , 


(2lj.) 


,2^n-2a 
d  b 


(J 


n/r  V)  =1 x__„  ^3+l^n-2a^2   ,  n 


-2a  1 


d"S 


;-?7n-2a 


n,.  n' 


d^ 


2^9 


d? 


±^— )^+(p^W,)  (25) 


if  2a  <  n  <  Ij.a  . 


It  may  now  be  assiriied,  witli  no  loss  of  generality,  that 
Pq^(5)  t   0  and  v^-f  {ii)   ^  0.   This  e::^sentially  determines  the 
value  of,  a,  and  limits  the  surface  forces  to  those  i-jhose 
orders  of  magnitude  o.rc  nob  less  than,  e^'^  =  h/R.   Thus  eqs. 
(22)  and  (2!).)  are  the  differential  equations  iiecessary  to 
determine  the  first  2a  non-vanishing  interior  stress 
coefficients. 
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To  surm'iarize  case  I-A,  ue  have  obtained  differential 
equations  for  the  first  2a  orders  and  stress  displaceraent 
relations  for  each  of  the  fir:rt  2a  non-vanishing  stress  coef- 
ficients.  All  these  results  are  deduced  vjithout  reference  to 
the  boundary  conditions  (10b)  along  the  edges.   VJe  shall  now 
explicitly  list  the  differential  equations  and  stress  displace- 
ment relations  of  the  tvo  lov/est  order  Interior  problems.   ?rom 
(22)  and  {2l\.)    employing  (20)  and  (21)  the  differential  equations 
of  the  "zeroth  order  interior  problem"  are: 

X  _  o  U   ,  ^,  dVJ   _  „  \ 


°^     dC" 


d? 


(26) 


— ?-3Sc p[ i,-+3VJ  +3v-^-J  =  3(Pt  -P, 


while   the    stress   dlsplaceraent   relations    (20),    (23)    and    (25) 
may  be  urltten   as, 


l._rc5U^^,.,.o      dV 


^  l-v"^    -^'^  dS 

a  2  o 

,-.0f,-    yi    _       1     r,.dU      -TiO      ,,d  VJ    ^1  j 

^  1-v^      ^^  d^"  --' 

a    /;.    y\  1  d  W    /,    ^2v  '  ■■^, 

^^  2(l-v")      d^-^ 


(27) 


,2a,,,    „v    _       1     r    1   d^VJ°y3   v   d^W°y2./l   d  •¥     .  ,,o  .„dU^  ^  ^t  .  ,    2a^   2a  v 
^  1-v^     ^     d£'+       "^     d^^-  "^     d£^  ^^  10    . 


(28) 


IS 


Equations  (26)  and  (2?)  are,  in  our  notation,  the  differential 
equations  and  stress  dicplacei.ient  relations  of  the  classical 
theory  of  thin  shells  due  to  Love.'^'     The  relations  (28) 
are  not  £;iven  in  the  classical  theory.   Si'iillarly  the  dif- 
ferential equations  rnd  the  stress  displacement  relations  of 
the  "first  order  interior  problera"  can  be  obtained  from  eqs, 
(20)- (25).   They  are  identical  vjith  those  given  in  (26j)-(28) 
provided  that,  1,  is  added  to  the  suDerscri^t  of  each  term. 

We  call  the  result,  eqs.  (29).   The  resultin.^^  stresses,  o''j 
a+1 


1        a+1    T  _2 
0.,  cr-^^  and  cr^ 


,  which  are  determined  from  the  solutions 
of  (29)  provide  the  first  interior  correction  to  the  zeroth 
order  interior  problem  consisting  of  tlie  classical  thin  shell 
theory  (26)  and  {2J),    and  the  additional  relations  (28). 

To  complete  the  formulation  of  the  interior  problems 
ap;:)ropriate  boundary  conditions  are  derived  in  Sections  li   and 
5  for  the  differential  equations  (26)  and  (29). 

VJe  continue  in  this  fashion  to  obtain  differential 
equations  and  stress  displacement  relations  for  the  higher 
order  interior  problems.   For  example,  for  the  2a+l  interior 
problem  we  obtain, 


3        2 
u3^(§,^)  =  u3^(^)^^  +  U^^(^)^-4-ij-3a(^)^+u3a(^j  ^ 

3. 


cr^i^.K) 


sj^(^);^^+sj"(^)^"+sj^(^)^+sj'^(?)  , 


l^{t,,K)  =   s^^(5)^^+s2^(c)^-  +  s^^(^)^H-s^^(C)  . 


(30) 
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The  expressions  for  the  coefficients  of  ^  ,  i  >  0,  in  (30), 
which  are  functions  of  W°(.?),  V/^^(?),  a^(^),  Tp^^'l.?),  Pq^(>^) 
and  ipji^)    will  not  be  explicitly  given. 

This  concnudes  the  discussion  of  Case  I-A.   Case  I-B  for 
vjhich,  2a-b  >  C,  is  analyzed  in  a  similar  raanncr.   Soi'ie  of  the 
results  are  given  in  the  appendix,   lie  observe  that  the  zoroth 
order  interior  stress  coefficients  can  be  detcrralned  only 
within  an  arbitrary  linear  function  of  E,, 

We  define  a  "reasonable  interior  shell  theory"  as  one 
whose  differential  equations,  boundary  conditions  and  stress 
displacement  relations  lead  to  uniquely  determined  interior 
stress  and  displacement  coefficients  of  every  order.   In  terms 
of  this  definition,  Case  I-B  yields  an  unreasonable  theory  and 
is  discarded.   Thus  we  have  shov;n  that,  viithin  the  limits  of 
applicability  of  the  expansion  procedures  presented  in  this 
paper,  there  is  only  one  scaling  of  the  coordinates  (16) 
vjhich  can  give  a  reasonable  theory. 


I|.,   Formulation  of  the  Boundary  Layer  Problem 

It  was  shoii/n  in  the  previous  section  that  the  boundary 
values  of,  o'"^  and  <t^  .s  are  explicit  functions  of  ^  [see  eqs, 
(20)  and  (23)]  while  the  applied  edge  stresses,  T(^|e)  and 
'ziK,;^),    are,  within  the  limits  of  equilibrium  and  continuity, 
arbitrary  functions  of  ^,   This  indicates  that,  in  general, 
the  interior  exoansions  cannot  satisfy  the  boundary  conditions 
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(7b)  along  the  edg:e.   Thus  if  the  series  (9)  represent  the 
solution  of  the  exact  theory,  tney  do  ao  in  a  region  awaj?- 
frora  the  edf^e ,    i.e.  in  the  "interior''  of  the  shell.   The 
expansions  therefore  cannot  hold  uniforraly  up  to  the  ed^e  and, 
in  factj  are  not  valid  at  the  ed^e.   The  solutions  of  the  exact 
theory  rrpidly  deviate  frora  the  expansions  (9)  upon  ap'proaching 
the  ed^e;  the  smaller  e  is  the  iiiore  rapid  the  deviation.   The 
rer'ion  in  which  the  rapid  deviation  occurs  is  called  the 
boundary  layer.   In  the  linit  as  e  — s*  0  the  "vidth",  5,  of 
this  boundary  layer  ->  0.   To  obtain  an  approximation  of  the 
rapidly  varying  stress  distribution  near  the  ed^e  we  employ 
the  teclinique  ''  ''  '"^  of  expanding  or  "stretchin;"^  '  the  variable 
t,    so  that  the  resultin,^  expansions  uniformly  renresent  the 
solution  up  to  the  edge.   If  we  introduce  the  "boundary  layer 
coordinate",  X,  defined  by, ' 

where  (2a)  i;;;  used,  then  in  the  lirait  as  e  — >  0  the  vjidth  of 
the  boundary  layer  which  vanis]:;ec  when  expressed  in  the  E, 
variable,  reraains  of  non-vanishing  length.   In  fact  any  fixed 
neighborhood  of  the  edge  in  the  t,   variable  corresponds  to  an 
arbitrarily  large  neighborhood  in  the  ni  variable  provided  s  is 
made  sufficiently  small. 

The  scaling  and  stretcliing  of  the  axial  coordinate  intro- 
duced in  (2a)  and  (31)  yields  a  sjrsteuatic  method  of  studying 


7  Here  we  are  examining  the  boundr.ry  layer  edjacent  to  the  edge 
^  =  x  =  0.   An  identical  procedure  applies  to  the  edge  x  ~  L. 
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boundary  layers  of  dirfcrent  xddths.   In  fact  if  we  define  the 
order  of  iaarniitude  of  5  as  the  order  of  raagnitude  of  the  dif- 
ference between  the  sea]  in ^f^-  pnd    stretching  lengths  of  the  aj:ial 
coordinate,  :c,  then  froi,:  (16)  and  (31)  it  follows  (for  Case  I-A) 
that, 


6  =  0[s/Rh  (1-v^)]  . 


To  investigate  the  boundary  layer  it  is  convenient  to 
consider  formulation  B  of  Section  ,^  consisting  of  the  equilib- 
rium equations  (/.).)j  the  co'aoatibility  equations  (6)  and  tlie 
boundary  conditions  (7).   The  coordinate  '<',  eq.  (31),  is  intro- 
duced into  these  equations  and  solutions  are  then  sought  in  a 
form  similar  to  (9).   These  solutions  must  satisfy  the  boundary 
conditions  along  the  edges  and  surfaces  of  the  cylinder.   In 
addition  they  must  "match'',  in  a  ;:iannsr  to  be  precisely  for"nulated, 
with  the  interior  solutions  ai;ay  from  the  edge.   Appropriate 
boundary  conditions  for  the  differential  equations  of  the  interior 
problems  are  also  obtained  from  this  boundary  layer  investigation. 

VJe  define  the  "boundary  layer  stresses",  fC^^^^e),  as, 


'C   (x,r)  -t      (x.r) 


(32; 


z^'^''       E    *         "zx 


and. 


_.A,~(^2,^;e)  -  f^  -1-  f^  +  f^  . 


It  5-S  assumed  thst  ecch  boundary  layer  stress  component  can  be 
represented  by  a  power  series  in  e  oT  the  foria: 


fi'i,K;^)-"  Z^io  ^^('^^^)e^  . 


(33) 


where   f^  =   0   if   i  <   0.      The   qucntities   f^{^,^)    are    called   the 

boundary   layer    r,tr"-ss    coefficients.      Sub.stitutinfr    (31),     (3?) 

and   expansions   of  the    Corra   (33)    iiito    ('4)    rnd    (6)    and    equating 

coefficients  of  like  povers  of  e,  v;e  find  from  the  coefficients 
n . 


or   e 


-i^n 

.1.    + 

,n 


Vjf 


n 


jf  -^Z     (f^  -  f^)  -  0  ,       (^ 


^^zx       '=^^x        — '       i 


c^^  .9^ 


(3i|-) 


^  X  1  +  V  ,,,2  ■^~  ,;^ 


\ 


?2^n 


■^    ^z        1+v      ^,^2 


5-    +    n'     ~#    -    2T"        (f^    -    f  M     :=    0    , 


c^^ 


r-72f,n        XT-'     r 


(35) 


"-zx   ^   1+v   d'C\^>L         ^        9K  ^      ^ZX 


Here,    the   binomial   expansions   of   3(^;)    and    3    (^)    are    einploved, 

r— ' 

>    IS  defined  in  Section  3  and 


z:  A 


1  — 
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Boundary  conditions  for  the  f-^{^,c,)    on  ^  =  +1  are  obtained  by 
as.suraing  that  p^(5)  and  Pj(5)  can  bo  expanded  in  Toylor  series 
of  the  form; 


CO   pV3  =  v-"   pi>;Jpi^ 


wViere    (31)    is   used   and, 

U     d-^'p^lo) 


1=01    ... 

-J*      di:-^ 


0  ,  j    <    0 


Employing   this  result    and    (10)    vie   find   from   (31)-(33)    the 
boundary   conditions    for    f    {^L,K,)  - 


f"  (0,^)  =??(«)   .        f2(o,?)  =?"(?;)   .  (36b) 

To   com.plete   the   fornulation  of   the  boundary  laysr  proble-.n 
consisting  of    (3lj-)-(36)    vje  reqtiire    that    the   boundat^y   layer   stress 
coefficients  match  the    interior    stress   coefficients    in   bhe   fol- 
lowing manner.      Assuiaing   that,  n:    (C,^),    are   regular   functions   of 
■5   at   £,   =   0,    xve   expand   them   in   a  Taylor   series    in  a   neighborhood 
of   the    edge, 

oH^,K)  -T.%0  'Y^^^^'  ' 

■where. 


1.  ^^'rHo,K) 


i  <  0  . 


It  follox^;^  fror:i  this  result,  (9)  £^nd  (31)  that  in  a  neichtorhood 
of  the  ed;;^e, 

a-(5,?K  z:T=o  nf^o  ^Yzyih'"-^^  =  E^.o  o'^^'i^^^^  .       (37a) 

Here, 

are   the    interior    sti^ess    coefficients   near    the   boundajry   c.s 
functions    of  ^l   £.nd   ^. 

From  the   definitions    (2ta)    and    (32)    it   fo],lox-js   that, 

o-(?,r,;e)   =  f('^,^;e)    . 

Employing  the   re^"^  re  sen  tat  ions    (33)    and    (37a)    vfe   tlien   find  the 
matching   conditions    (or  the    as^n'ivototic   form)   for    the    stress 
coefficients   for   large    ^  .      V'e   virite   these    as, 

lim       f"(/?.,^)   =    c-^^i^L^K)    ,  n  =  0,1,...,  (38) 

Before    exaraining  the   system    (3'4)-(36)    snd.    (38)    in  detail 
a  result    is   stated,    wiLhout   Di'oof,    concea^ning  the   boundary  value 
problem  for   a  function,    ^{^,^),    defined   on  the    semi-infinite 
strip,    -^  >   0,     1^1    <   1,    and   satisfying   tlie  biharraonic   equation. 
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'■  ■  =  2i^  +  o  J^_  ^  ?JH  = 


-Z^     ■  ..y..^^- 


and  the  boundary   conditions, ' 


^ 


^yy{0,-C,)    =  h^(^)    ,      4>-^.(0,^)    =  hp(^)    ,      c|v^(x,+l)    =   4/^^(::,J:l)    =0,! 

}(39) 

liri      (4     (^6,^)  ,     ^^.A\,r^))  =  0  .  ) 

It   is   easily  shov;n  that    if   4vv»    4/,«--   ^'<^^  4v   a-^''-    uniforrrj.y   con- 
tinuous  functions    of  ^   as  /^   — 5»  cd  ,    then  in  order   that   4}    4  a 
and  4v   sre    single   valued   functions, 

1  1  1 

-1  -1  -1 

These  results  vdll  frequentl"  be  referred  to  as  the  "integral 
relations". 


5.   Analysis  of  the  Boundary  Layer  Problera 

i/Je  consider  Case  I-A,  b  =  2a.   Iquations  (3l|)-(36)  yield, 
for  n  <  Za, 

.-.^n   ,-^^.n  -,^:,n    ^^:,n 


'   Subscripts  on  4('^>^)  indicate  "partial  differentiation,  e,,^,,, 

4   =-^fi 
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,2  A  n 


(^4: 


V72..n  _   ,.  V72f>"     +  _J-_  ^  -/^    =   o 

f^,^J''l,+l)   =  f^(>?.,+l)  =  0   , 

To   apply   the   matching   conditinns    (3?;)    we  r-ust   consider   tv/o 
separate    esses    for   tb.c  ran-^-e  n  <   .?a;      (i)      0  <  n  <   a, 
(il)      a  <   n  <   2a.      Note    thr.t    if   the   external  forces   en   the 
surfaces    are   of   order    e"  then  a   =  1   end    (i)    md    (ii)    becorae, 
n  =   0   and  n  =  1,    respectivel-^-^.      '.'e   find   usino:   the    interior 
stress    coefficients    (17),    (20),     (21),    (23)    and    (37b)    that   the 
matching  conditions    (38)    are  wi-itten  as, 

[S;^(0)+s5(0)4'    ,    0    ,    s2(0)+3^:^(())^    ,    0]    for    (i 


(^2) 


X 


dS-2"^(0) 


lira     [f^,f^,fS,f^^J  =/'[s^(o)+sr;(o)^+~'V^ ?^^  ,  0  ,  s2(o)+3^(o)^  y(i;)4) 


i^-SoOO 


ds2~^(o)     ds-r^'lo)  .  dsJ"^(o)        p 

for    (ii)    , 

where,    (22)    is   also   ucod   and    an   ohvicus   notation   is   eraployed. 
Solutions   rnsv  be   obtained   to    ( ','1 )- (!;.'|)    by   defining  stress 

functions   Y^(^,?)    on   the    serai-infinite    strip,    ^i  >    0,     |.>^|    <   1, 

t 

such  that, ' 


'      Subscripts    on  y(^,^)    indicate   partial  differentiation. 
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and  r  V\S) 

fg  =  v(fx  ■"  ^S^  "  ^"^^^^Y  ,     0  <  n  <  2a  ^1 

By  direct  substitution,  xje  see  thct  ('i.5)  are  formal  solutions  of 
(/|1  )  and  (lt-2)  provic'ed  that  y  (%^)  -"'^-  biharrionic  functions. 
To  a'p:ply  the  inte,;:ral  i^elrtions  (li.O)  x'c  define  the  reduced 
stress  function,  cj/  [^i,c,),   and  the  reduced  boundary  layer  stress 
coefficients,  P^'(^.,^j,  as, 

T,n    in  _  n    t  .    .^n    T-,n  _  in  _  n    ,  .    ^^n   1 

^       ^    'i— >C0  O-'8!>00   ^ 

It  then  follows  from  {bU) ,    ik^)    and  ([l3)  that  the  ^^    satisfy, 

VH""  =  0  ,  0  <  n  <  2a  ,    (I.?) 

and  assume  the  boundary  values  ('i.O)  i;here, 


h^  iK)    =   ^''(^)  -  SiJCO  -  -;(0)^  ,   if  0  <  n  <  2a  , 

if  0  <  n  <  a   ,  y     {I4.Q) 


[ -  ^J(^)    ,  if  0  <  n  <  a   ,   V 

-  'c?(^)+-o  — "^v (1-r)  .   if  a  <  n  <  2a  .   / 


h^ia 


Substitutlnp-  {I4.Q)    into  the  integral  relations  ([j.0)  vje  find  that. 
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Sj(0)  =  I  I  a^{K)c\K    ,  I  (^-!-9a) 

"^-  ^  /  0  <  n  <  2a 

sn(o)  =  .   1   dflAoj,  ^  3   f  ^.f^-(^)d^  ,  (i+9b) 


1  — », 

j  r'^{^)dK   =  0  ,  0  <  n  <  a  ,   "j 

-1  m9c) 

It  follows  from  the  first  of  (I|.9c)  that  if  an  edge  shear 
force  is  applied  of  order  of  magnitude  lesr,  than  s  ,  (see  (10b)), 
then  it  must  be  self-equilibrating.   In  the  follov;ing  we  restrict 
ourselves  to  those  ^^(^)  and  j^(^)  that  are  not  self- equilibrating, 
Then  without  loss  of  generality,  vie   assume  that,  o*  (^)  ^  0, 
t-^(?)  ^  0,  ~c^^^  =   0  'J-^  n  <  a.   Equations  ()|9a,b)  then  yield 
two  of  the  boundary  conditions  for  the  first,  2a,  interior 
problems,  while  the  second  of  (';-9c)  yields  the  remaining 
condition  for  the  first,  a,  interior  problems.   The  tlir-ee 
boundary  conditions  for  the  differential  equations  of  the 
zeroth  order  interior  problem  are,  from  (i4-9), 


S^(0)  =  -^    J  c>    (^)d^  , p~^   =  -  -"-^ '-     )  <<:•  (^)d4  , 

-1  "^^  -1 

dV(g)  _  W)  |^(^),^  . 


(5o: 


d^  _i 
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Equations  (5o)  arc,  in  terras  of  the  variables  of  this  paper, 
the  boanclary  conditions  usuall.v  assm-aecl  in  the  claosical  thin 
shell  theory  of  Love.  *     If  a  >  1  then  the  boundary  conditions 
for  the  first  order  interior  nroblem  raajr  also  be  obtained  from 
(1^9).   They  are  identical  xxTith  those  given  in  (5o)  provided 
that,  1,  is  added  to  the  superscript  of  each  terra.  Vie    call 
these  boundary  conditions,  eqs.  (5l).   If  a  =  1,  ([|.9a,b)  vill 
give  the  first  tuo  boundary  conditions  \ov   the  first  order 
interior  problem.   The  renainin^a'  condition  is  obtained  from  a 
similar  analysis  of  the  boundarjr  layer  equations  (l|J)-([i.3)  and 
the  matcl'iins  conditions  for  n  =  2a.   These  calculations  are 
quite  lengthy  end  are  not  exliibited.   riov:ever  they  shov;  that 
this  boundary  cond5,tion  coincid^.s  with  the  third  of  (5l).   Thus 
(^1)  are  the  boundary  conditions  fox^  the  differential  equations 
(29)  of  the  first  order  Interior  problem  for  all  values  of,  a. 

It  is  Important  to  observe  that  the  prooer  boundary 
conditions  for  the  interior  problems  are  obtained  by  investigating 
the  boundary  layer.   In  addition,  this  investigation  yields  a 
sequence  of  boundary  value  oroble-.,is,  con":  1st  in f?  of  (l;.?),  (39) 
and  ([|-3)  to  systematically  calculate  the  stresses  in  the  neigh- 
borhood of  the  boundary.   For  convenience  ve  list  here  th.e 
first  two  "boundary  layer  problems'' j   the  zeroth  order  boundary 
layer  problem. 
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V 'cr  =  0 , 


4;°    (^^+1)    =   4^/,(^?,+l)    =   0    ,      llr,i    [(i)°^('?,.:),4^^(^^,^)]    =   0    ,  ^  (52) 


1 


z^— >CO 

1 


-1-1 

the  first   order   boundary  layer  problera. 


V^^i  =  0  , 

U^  '  /([  — >co     ^  =^ 

1  1 

4jv(o,^)  =<F^(^)-iJ~^(^)d^-|^  j<^~-^(^)d^  , 


-1 


4j^(o,^) 


(  -<Vl(?)+^(l-?^)  j   t;];{^)d^    ,      if    a   -   1    , 

-1 


(53) 


,      if   a  >   1    .  / 


V. 


6.   A  Summary  of  the  Shell  Theory 

By  the  terra,  ''thick  shell  theory''  we  imply  a  prescription 
for  iiiodifylng  the  results  of  thin  shell  theory  to  obtain  a  more 
accurate  a'opro:ciraation  of  the  stresses  In   a  shell  considered  as 
a  tliree  diraensional  elastic  body.   Vie  expect  a  thick  shell 
theory  to  su-oply  the  doninant  corrections  to  the  results  of 
thin  shell  theory  over  the  entire  cylindrical  region. 

In  the  preceedlng  sections  \:e   described  a  netliod  of 
obtaining  fro^a  the  exact  theory  systeras  of  differential 
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equations  and  ntress  displrceraent  relation.^  vhose  solutions 
systematically  aoproicir.iate  the  stresccs  In  the  cylinder.   Two 
tyoes  of  e::pansions  in  terras  of  e,  are  employed:   The  first  of 
these,  the  Interior  expansion,  provides  approximations  of  the 
stresses  in  the  interior  of  the  shell;  while  the  second,  the 
boundary  layer  expansion,  yields  a  description  of  the  stresses 
near  the  ed^e.   Appropriate  boundary  conditions  for  the  dif- 
ferential equations  that  result  from  the  interior  expansions 
are  deduced  from  properties  of  the  boundary  layer  expansions. 
V'e  now  sui'imarize  the  results  of  Sections  3,  If.  and  5  to 

define  a,  ''thick  shell  theory  of  order  N".   Let  the  stresses  of 

th 
the  H   order  thick  shell  theory,  represented  by  the  ce^O-ric 

( F ) 
symbol,  '^   (?,?,£),  be  defined  as, 

where  cr      and  F   are,  respectively,  the  interior  and  reduced 
boundary  layer  stress  coefficients.   Thus  the  thick  shell  theory 
of  order  zero  for  which  -'  '   =  c     +  F  ,  is  composed  of  the  zeroth 
order  interior  problem,  eqs,  (26)-(28)  and  (5o),  and  the  zeroth 
order  boundary  layer  problem,  eqs,  ( L[.6 )  and  (52).   Therefore 
this  thick  shell  theory  modifies  the  classical  thin  shell  theory, 
eqs.  (26),  (27)  and  (5o)  by  supplyin,^:  a  first  approximation, 
valid  in  the  interior,  of  tVie  transverse  shear  and  normal 
stresses  (28),  and  a  cori-ection  which  results  froii  (''6)  and  (52), 
of  the  stresses  in  the  nei'rhborhood  of  the  edge.   In  a  similar 
manner  the  thick  shell  tleory  of  order  one  for  v-'hich 
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cr^  '  =  T  +  F   +(o~  +7  )e,  ic  obtained  by  corablning  the  results 
of  the  zeroth  order  thick  sliell  theory  v.'ith  the  contributions 
of  the  first  order  interior  problcra  (29)  and  (5l)  and  the  first 
order  boundary  layer  problem  ('|6)  and  (53).   VJe  can  proceed 

in  this  fashion  to  obtain  thick  shell  theories  of  any  order. 

15 
Johnson  and  Reissner   proposed  corrections  to  the  results 

of  thin  shell  theory  uhich  i-ere  obtained  by  an  ''interior 

exToansion''  -procedure  th^t  is  similar  to,  but  less  ri;enei'al  than, 

ours.   If  v;e    specialize  cur  results  to  the  particular  ;oroblein 

treated  in  ref,  l5  (for  the  icotropic  case),  it  can  be  shovn 

that  their  first  correction  coincides  ^rlth  th?t  furnished  by 

our  second  order  interior  probler,!  (30)  and  is  therefore  0(e  ). 

Since  the  boundary  layer  is  not  analyzed  in  ref.  1$,    the  lower 

order  corrections  resulting  frova  the  zeroth  and  first  order 

boundary  layer  problems  and  the  fir-it  order  interior  problem 

are  omitted. 
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/roDenc'lx 


Case    I-B,    ?a   -  b  >   0;      We    can   show  from   (ir%}--(13)    that 


w    (4,^)    -  W    (s)    =   c^'^    +   c^    , 


.A+n, 


u        (5,^)   =   -   c^'^,      if   n  <   2a--b,    or   X+n  <   a      , 
v.^{t,,a    =   -  ^^^1  ^   +  U^(a    ,      if  n  =   a    , 


vjhere, 


t;;(?)  =  . 


1     6^\r^'^U) 


Tn'-'  '  r   \      _ 


1-v^  d^^ 


Bj'-  proceeding   in  a  fashion   si-iilej-   to   that   'ased   in  Case    I-A,    i 
can  deduce,    assuming  p^"^  =   0   if  n  <   2?.,    that. 


-/e 


cr       =0 
zx 

^n     _    „ 


if  n  <   X   , 
if  n  <   2\   , 


'-'zx        2  "d?  ^^"^    ''    ' 
^z  2       ,72^^      T^^o    ^Pl      ' 


X     _  1        rd     U       ,„^0l      _     „ 


J  2;sO  , 


.    /S    :■  ii 
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where  c^  and  Cp  are  arbitrary  constants.   The  above  equations 
forra  the  stress  displacement  relationG  and  the  differential 
equations  for  tiie  zeroth  order  interior  protlera  of  Gace  1-3, 
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